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Abstract In the framework of the three-party constrained voter model, where 
voters of two radical parties {A and B) interact with "centrists" (C and C^), we 
study the competition between a persuasive majority and a committed minority. 
In this model, A's and B's are incompatible voters that can convince centrists or 
be swayed by them. Here, radical voters are more persuasive than centrists, whose 
sub-population consists of susceptible agents C and a fraction ^ of centrist zealots 
C^. Whereas C's may adopt the opinions A and B with respective rates 1 + 5a 
and I + ds (with 5a > ds > 0), are committed individuals that always remain 
centrists. Furthermore, A and B voters can become (susceptible) centrists C with 
a rate 1 . The resulting competition between commitment and persuasion is stud- 
ied in the mean field limit and for a finite population on a complete graph. At 
mean field level, there is a continuous transition from a coexistence phase when 
^ < Ac = Sa/ {I + Sa) to a phase where centrism prevails when > Ac . In a finite 
population of size A'^, demographic fluctuations lead to centrism consensus and 
the dynamics is characterized by the mean consensus time T. Because of the com- 
petition between commitment and persuasion, here consensus is reached much 
slower (^ < Ac) or faster (i^ > Ac) than in the absence of zealots (when T ^ A^). In 
fact, when ^ < Ac and there is an initial minority of centrists, the mean consensus 
time grows as T - TV" ^ /^e^'', with TV > 1 and y = 5a - C ( 1 + In ( 5a / C) ) + ^ ( 5| ) . 
The dynamics is thus characterized by a metastable state where the most persua- 
sive voters and centrists coexist when 5a > 5b, whereas all species coexist when 
5a = Sb- When ^ > Ac and the initial density of centrists is low, one finds T ^ InA^ 
(when 1). Our analytical findings are corroborated by stochastic simulations. 



Mauro Mobilia 

Department of Applied Mathematics, Scliool of Mathematics, University of Leeds - Leeds LS2 
9JT, U.K.; 

Tel.: +44-113-3431591 
Fax: -1-44-113-3435090 
E-mail: M.Mobilia@leeds.ac.uk 



2 



Keywords Opinion dynamics • Voter model • Zealots • Consensus time • 
Metastability 

PACS 02.50.-r • 89.75.Fb • 87.23.Kg 

Mathematics Subject Classification (2000) 60J28 91D10 91A22 
1 Introduction 

Opinion dynamics aims at understanding how cultural change evolves UHJllll. 
This issue is closely related to problems arising in various disciplines in the life 
and behavioral sciences f4\, and the recent years have witnessed a growing activity 
in applying the tools of statistical physics to their description |3,5|. In models of 
opinion dynamics, each agent is typically represented by its "opinion" state (often 
labeled as a spin variable |6|) that is updated in response to the opinion of a lo- 
cal neighborhood, a mechanism of cultural dynamics documented in sociological 
studies, see e.g. |7,8|. Basic issues in opinion dynamics concern the time neces- 
sary for a consensus to be attained and the period during which cultural diversity 
is maintained ll3]|3. These questions are usually addressed by considering simple 
and insightful models, like the influential two-state voter model (VM), in which 
local opinions are copied and spread by imitation, see Sec. 2 below. It has recently 
been proposed that the evolution of cultural diversity would be more realistically 
described by models where the seek for consensus is limited by some form of in- 
compatibility |[T]|2l. This idea has been investigated in an analytically amenable 
three-party constrained voter model (3CVM) IWlQlTH I (see also Ref. 1 12|). 

In this work, we generalize the 3CVM to study the competition between a 
persuasive majority and a committed minority. In our model, as in Refs. Il9lll0l 
[TT|, the voters of two radical parties (A and B) are incompatible and do not interact 
among them, but they interact with a third species ("centrists", C and C^). Here, 
radical A and B voters are considered to be more persuasive than centrists, whose 
sub-population comprises susceptible agents (C) and a small fraction ^ of "centrist 
zealots" (C^). The latter are committed individuals that always remain centrists, 
while a susceptible centrist adopts the opinion A with rate 1 + 5^ when it interacts 
with an A-voter, and the opinion B with rate 1 + 5b when it interacts with a B-voter. 
Furthermore, each A and B individuals can become a (susceptible) centrist with 
rate 1 . The parameters ^4 > 5^ > hence represent the persuasion biases toward A 
and B, i.e. the respective persuasion strength of these radical voters. The dynamics 
is thus characterized by a competition between the persuasion strength of A and 
B voters opposed by the resistance of centrist zealots (C^). The questions that 
we address with this model are the following: (i) When and how does a committed 
minority prevail against a persuasive majority ? ( ii) How does the mean consensus 
time vary with the fraction of zealots and the persuasion biases? 

These questions are answered by studying the generalized 3CVM's properties 
in the mean field limit and for a finite population of size A' on a complete graph. 
In the absence of fluctuations, we find that more than one opinion coexists when 
^ is below a critical threshold Ac — 5a/(1 + 8a), while centrism prevails when 
^ > Ac. The dynamics is markedly different when the size of the population is 
finite and the demographic fluctuations drive the system into centrism consensus 
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in a time that depends non-trivially on 5a , 5b , and A'. The mean consensus time 
in a large population with an initial minority of centrists is here shown to grow 
dramatically as T 7V-i/2gA'r(^A.?) when < A^., and as T InA'^ when ^ > Ac. 
These results, derived using the Fokker-Planck equation and the WKB approach, 
are in stark constrast with the mean consensus time obtained in the absence of 
zealots centrists (when T ^ N IIIOIIIU ). Our findings, corroborated by stochastic 
simulations, thus demonstrate that zealot centrists either significantly slow down 
< Ac) or speed up > Ac) the approach to consensus. 

While some aspects of the influence of committed agents in opinion dynamics 
have been considered in the literature, e.g. in Refs. |[T3][T4l[T5 II 1 6 IfTTl (see Sec. 2), 
this work differs from earlier studies in various respects: (i) We here consider a 
three-party (four-state) model whereas most of previous works, like Refs. 1131 
[T4l[T5l[T6l[T7]| . focused on two-state systems, (ii) Here, we study the competition 
between the degree of commitment in the population and the persuasion strength 
of two types of voters (three independent parameters), while no persuasion bias 
was considered in 1151 16, 171. (iii) Furthermore, whereas most of the previous 
results were obtained in the mean field limit |16| or by means of numerical simu- 
lations 1 17], the exponent 7(5^,^) is here obtained analytically. 

This paper is organized as follows: Next Section is dedicated to a brief review 
of the voter model and some of its variants, while the 3CVM with centrist zealots 
is introduced in Section 3. The mean field dynamics of such a model is discussed in 
Section 4. The section 5 is dedicated to the stochastic formulation of the dynamics 
in a finite population. The mean consensus time in the case of identical persuasion 
biases is studied in Section 6. In particular, the case of the long-lived coexistence 
of opinions and metastability is studied in terms of the Fokker-Planck equation 
(Sec. 6.1.1) and with the WKB method (Sec. 6.1.2). The mean consensus time in 
the general case of asymmetric persuasion biases is analyzed in Section 7. In the 
final section, our findings are summarized and our conclusions presented. 



2 A brief review of the voter model and some of its variants: zealotry, 
commitment and incompatibility 

As a contribution to a special issue on the "applications of statistical mechanics to 
social phenomena ", this work is dedicated to the study of the competition between 
commitment and persuasion in the three-party constrained voter model. Before 
presenting our new findings, it is useful to begin by providing some background 
material in the form of a brief review of some basic models and results. 

Since the pioneering works by Schelling and Granovetter fTT], the impor- 
tance of relying on individual-based models to relate "micro-level to macro-level 
interactions" |[8l in social dynamics has been recognized. As an inspiring example 
of this type of modeling approach, Schelling showed how homophily, that is the 
tendency of an individual to being bond with neighbors sharing similar character- 
istics, is important to understand the formation of social segregation [7 |. 
In this context, the voter model (VM) [61 is one of the simplest and most paradig- 
matic individual-based models of opinion dynamics. The VM describes how a 
socially interacting population of individuals possessing a discrete set of states 
("opinions") evolves toward consensus as the result of interactions between neigh- 
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bors and local fluctuations of the population composition. The VM is closely re- 
lated to the Isingmodel of statistical physics with the Glauber kinetics at zero tem- 
perature [18 19 [3, and also to the Moran model |20| commonly used to describe 
evolutionary dynamics in the life and behavioral sciences |4 0. In the classic two- 
state VM, each node of a graph is occupied by a "voter" in one of the two possible 
state, e.g. denoted A or B (or +/— , t / i> or "leftist/rightist"). These voters can 
be interpreted has having no self confidence and changing their state by simply 
adopting the opinion of a neighbor. Hence, in its simplest form the VM evolves 
according to the following rules: 

- Pick a random voter. 

- The selected voter adopts the opinion of one of its random neighbor. 

- These steps are repeated until a consensus (where all individuals are of the 
same opinion, either all leftists or all rightists) is necessarily reached. 

In the language of statistical physics, these rules and the VM define a Markov 
chain with absorbing states that is aptly described by a master equation |22|. One 
of the most appealing features of the VM is its mathematical tractability. Indeed, it 
is one of the rare models of non-equilibrium statistical physics to be exactly solv- 
able in any dimensions on regular lattices, and many of its properties are known 
analytically, see e.g. ll3l [T9ll22j . The most significant features of the VM are (i) 
the probability to attain a specific consensus state, (ii) the mean time to reach a 
consensus, and the two-point correlation function ?#(r,f). 

i. The consensus probability Ea ix) (also called "exit probability", or "fixation 
probability" by analogy with the evolutionary dynamics literature [4]) is the 
probability that a finite system consisting of an initial density x of A-voters 
eventually reaches the consensus state where all voters are in the state A. As 
the average opinion, also called "magnetization", is conserved by the VM on 
any degree-regular graphs (like lattices and complete graphs), and since the 
VM always reaches consensus, one has Ea{x) = jc on any regular graphs. 

ii. On regular graphs, it is known that the mean consensus time T (also called 
"mean exit time" or "mean fixation time") in a population of A' voters scales as 
X m one dimension {d = \) and as T NlnN in two dimensions {d = 2), 
whereas in dimensions d > 2 one has T ^ N (complete graphs correspond to 
d = oo). This implies that d = 2 is the upper critical dimension of the VM. 
Moreover, it has recently been shown that the mean consensus time scales 
sublinearly with A' on heterogeneous graphs characterized by nodes of high- 
degree, i.e. T ^ Af", where < a < 1 depends non-trivially on the network's 
degree-distribution ||2T1| . 

iii. Significant insight into the spatial distribution of opinions is provided by the 
two-point correlation function 'i^{r,t). This quantity informs on the probability 
that two voters separated by a distance r are in the same state at time t. On lat- 
tices of dimension d, ^#(r, f — S- oo) decays asymptotically as r'^~^ in high dimen- 



These models are equivalent on one-dimensional lattices, but not in higher dimensions 
where the Ising-Glauber model follows a majority rule, whereas the VM evolves according to a 
a proportional rule. As a consequence, while there is a surface tension in the two-dimensional 
Ising-Glauber model this is not the case for the VM, see |,3>19 l and references therein. 

^ The voter and Moran models are equivalent on regular graphs but lead to markedly different 
dynamics on heterogeneous networks, see e.g. \2\\. 
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sions (d > 2). In low dimensions the spatial organization of the VM is charac- 
terized by coarsening, i.e. the slow formation of growing domains of a single 
opinion. As a consequence, the two-point correlation function asymptotically 
approaches the value 1 in low dimensions. More precisely, when f — > oo and 
< r < V?, one finds 1 - ^(r,t) ^ rjsftm d = 1 and 1 - ^{r,t) ^ Inr/lnf in 

J = 2 1221. 

While the VM has proved to be insightful and influential, it relies on oversim- 
plified assumptions like the total lack of self-confidence of all voters and the un- 
avoidable formation of a consensus. In fact, since Granovetter's seminal work on 
"threshold models" |8|, the influence of the population's heterogeneous response 
to stimuli on social dynamics is well established. In the opinion dynamics context, 
as a step toward the investigation of a simple model describing the dynamics of 
a heterogeneous population with different levels of confidence, we have proposed 
to investigate the two-state voter model in the presence of "zealots" 1 13 , 14 15 1. 
The zealotry is implemented by assuming that a small fraction of the population 
consists of committed individuals (called zealots) that either favor a specific state 
{A or B) as in Refs. II13II14I . or adopt an opinion that never changes as in Ref. lilSI . 
As a consequence of zealotry, the magnetization is not conserved in the VM with 
zealots, and the main features of these models are the following: 

- In Ref. {T3\ it was shown that on low-dimensional lattices (when d <2) a 
single zealot with a favored opinion (i.e. a voter with a bias toward one opin- 
ion, but whose state is not fixed) imposes the consensus in its favored state to 
an infinite group of voters, whereas in higher dimension the magnetization is 
non-uniform and decays with the distance from the zealot. Furthermore, the 
unanimity state imposed by the zealot is approached algebraically in one di- 
mension, as f~^/^, and much slower in two dimensions (^ 1/lnf)- It has also 
been shown that on a finite lattice of size N = (where L is the lattice lin- 
ear size), the time necessary for the zealot's opinion to spread is is d = 1, 
NlnN ind = 2, and A' in higher dimensions. 

- Ref. 11411 was dedicated to the study of the VM in the presence of a small group 
of n zealots, each with a favored opinion and a specific bias toward its preferred 
state. The competition between zealots of different types was shown to prevent 
the formation of consensus and to lead to a non-trivial fluctuating steady state 
whose properties were investigated by exploiting a formal analogy with the 
electrostatic potential generated by n classical point charges. The approach 
toward the steady state was shown to be algebraic (^ f~'/^) and logarithmic 

1/lnf) in one and two dimensions, respectively, and to be characterized by 
the formation of growing domains when the density of zealots is very small, 
i.e. n/N ^ 1. In this case, the size of the single-opinion domains was found to 
grow with the system size but to never span entirely the system. 

- The voter models with a finite (but small) fraction n /N of zealots as individ- 
uals that never change opinion has been studied in Ref. ifTSl . The study was 
carried out on complete graphs and on low-dimensional lattices with randomly 
distributed zealots (of both types). It was thus shown that in all cases a small 
fraction of zealots is effective in maintaining a reactive state characterized by 
a Gaussian magnetization distribution with a width that decays as n^'/^. 
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The influence of committed agents in opinion dynamics has also been considered 
in other two-opiruon dynamics models, see e.g. Refs. II16II17II24II25I : 

- The case of a committed minority of "inflexibles" in a two-state majority-rule 
model has been considered at mean field level in Ref. |FT6|, where it was shown 
that an equal densities of inflexibles of each type prevent consensus from being 
achieved. 

- The authors of Ref. 11711 studied how the mean consensus time varies with 
the fraction of committed individuals in a two-opinion variant of the naming 
game |26l called the binary agreement model. In the model of Ref. 1 17|, the 
population initially consists of a majority of individuals of one opinion and 
a fraction p of individuals of the other species. As in Ref. [15], the latter are 
committed individuals that never change state and impose their consensus after 
a mean time showed to grow exponentially or logarithmically with the popu- 
lation size, depending on whether p is below or above a critical threshold 1171 . 

- The authors of L24 1 considered the confident voter model, in which each voter 
has two levels of commitment ("confident" and "unsure") for each of the two 
possible opinions. The confident VM has been studied in the mean field limit, 
where the consensus time is T In A', whereas in one and two dimensions the 
mean consensus time scales respectively as T A^^ and T ^ N^/'^ 1241 . 

- Furthermore, the heterogeneous voter model in which each agent has its own 
intrinsic rate to change state has been considered in Ref. |25|, where it was 
found that the time until consensus is reached is much longer in the heteroge- 
neous VM than in the classic voter model. 

In addition to voter-like models with committed individuals, there are vari- 
ous prototypical opinion dynamics models, whose outcome is characterized by 
a lack of consensus. Influential models of this type are the multiple -state Axel- 
rod model [l] and the bounded compromise model |2|. The key feature prevent- 
ing consensus in these models is a form of incompatibility: when the opinions 
of two agents are too different, they are deemed to be incompatible and there is 
no interaction between such individuals. This can lead to cultural fragmentation 
and a frozen stationary state where a mixture of incompatible states coexist. In its 
essence, the feature of incompatibility is captured by the so-called three-party con- 
strained voter model (3CVM) that can be regarded as a discrete three-state version 
of the bounded compromise model ll9l[lOI. In the 3CVM, there is no interaction 
between species A ("leftists") and B ("rightists") that are are incompatible. How- 
ever, A- and B-voters interact with a third species C, and thus compete to impose 
their own consensus according to the following rules: 

- A centrist C can become an A-voter or a B-voter with rate 1 + 5' (and < 
|5'|<1); 

- A and B voters can become centrists C with a rate 1 ; 

- A and B voters do not interact: AB — > AB, 

where 5' is a bias favoring either A and B (when 5' > 0), or C (when d' < 0). 
The 3CVM admits three absorbing fixed points, corresponding to the consensus 
with A,B and C, and a "polarization line" along which a frozen mixture of non- 
interacting A's and B's coexist. This 3CVM has been solved analytically on a 
complete graph and the probabilities of reaching each consensus state and the 
polarization line were determined, both when 5' = 11 Oj and in the case d' 7^ 
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ifTll . Within an approach based on the Fokker-Planck equation |22, it was also 
shown that the mean exit times in the 3CVM scale linearly with the population size 
when the intensity of the bias is weak (i.e. when N\5'\ <^ 1 and A'^ ^ 1) ifTOl 
[Fl |. The 3CVM has also investigated on regular lattices and found to exhibit slow 
kinetics in low dimensions |9| (see also 112l ). 

Aiming to study how the maintenance of cultural diversity is affected by the 
competition between persuasion and commitment, in this work we consider a 
three-party (four-state) opinion dynamics model that combines the main features 
of both the 3CVM and VM with zealots. The resulting model, whose detailed 
specification is given in the next section, is here used to study the competition 
between two radical opinions (A and B) that are incompatible and the voters of a 
third party that are less persuasive but on average more committed than A's and 



3 The three-party constrained voter model with centrist zealots 

In this work, we introduce and study a three-party constrained voter model with 
centrist zealots on a complete graph (used for its tractability and because it is ar- 
guably a natural first choice in the context of social dynamics). Like the 3CVM |[9] 
[TOlim . this model is characterized by a population of individuals, j are of 
species A, k of type B and i are "susceptible" centrists (species C). Moreover, 
the population also comprises committed centrists (zealots), denoted C^, that 
never change opinion, as in II15II . Hence, one has A^ = j + + £ + with a fixed 
fraction ^ = /N of zealots C^. In the language of the voter model (see Sec. 2), 
the individuals of species A ("leftists") and B ("rightists") are voters of two radi- 
cal parties, while supporters of the third party ("centrism") are either susceptible 
centrists (species C) or centrist zealots (type C^). As in the 3CVM, radical opin- 
ions are incompatible and there are no interactions between A's and S's, but these 
voters interact with a centrist neighbor according to the following evolutionary 
rules: 

1 . A susceptible centrist C can become an A-voter with rate 1 + 5a and a 5-voter 
with rate 1 + 5b; 

2. A- and B-voters can become susceptible centrists with a rate 1 ; 

3. A- and B-voters do not interact: AB — )■ AB\ 

4. Centrist zealots always remain in this state. 

In the spirit of the VM and 3CVM, the dynamics of the system is thus implemented 
as follows: 

(i) At each time-step a pair of ("neighboring") voters is randomly picked; 

(ii) if an A/B-centrist pair is picked, the population composition changes accord- 
ing to the following schematic moves: 



S's. 



AC ^CC 



AC ^AA 



BC ^CC 



BC ^BB 



with rate 1+8^ 

with rate 1 + 5g 

and AC^ — > CC^ with rate 1 

and BC^^CCj^ withratel; 



(1) 
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(iii) if the randomly picked pair of neighbors consists of two radical voters {AA,BB 
or AS) or two centrists {CC,CCq or C^C^ ), the composition of the population 
does not change; 

(iv) these steps are iterated until centrism consensus is necessarily reached (see 
below). 

In this work, we assume the existence of persuasion biases (or persuasion 
strengths) 5a>Q and 5g > toward opinions A and B, respectively. The quantities 
^A/B measure the biases toward opinions A/B (against centrism) and thus reflect 
that A and B voters are more persuasive than centrists. Here, without loss of gener- 
ality, we assume that 5a > 5b > 0. The limiting case = 5* = = corresponds 
to the 3CVM without zealots and no bias that was studied in Ref. [ 10 1, while the 
3CVM with a symmetric bias but without zealots of Ref. OIU is recovered when 

From the rules ([l]), we expect a subtle competition between A and B voters, 
favored by their persuasion biases and seeking to impose their consensus, and 
centrist zealots that resist the spread of A's and S's and strive to promote centrism. 
We are particularly interested in such a competition when centrists are initially 
in the minority and the population is finite. In this situation, we will study the 
circumstances under which a committed minority can prevail against a persuasive 
majority. Another question of great interest concerns how the competition between 
persuasion and commitment leads to a long-lived (metastable) coexistence state 
along with the maintenance of a form of cultural diversity. These issues are studied 
in the mean field limit (Sec. 4) and in finite populations (Sees. 5-7). 



4 Mean field analysis 

At mean field level, one assumes a population of infinite size (A'^ — ?■ °°) and ignore 
any random fluctuations. Hence, the densities a = j/N,b = k/N and c = l/N, 
and the fraction C = of (centrist) zealots, are treated as continuous variables 
obeying the following rate equations (REs): 



ja{t) = a{t)[8Ac{t) - Q 

jb{t)=b{t)[8Bc{t)-Q. (2) 



The total population size being conserved, one has c(t) = I — ^ — +b{t)). 
Substituting this expression into Q yields two coupled equations whose proper- 
ties are discussed below, first in the case of asymmetric biases and then when the 
biases 8a and 5g are identical. In fact, we readily notice from (|2]l that in the ab- 
sence of bias {8a = 8b = 0), the system quickly evolves toward a population com- 
prising only centrists, whereas an interesting situation arises when 8a> 8b > 0. 
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Fig. 1 Schematic phase portrait of the mean field dynamics {2}. (a) Asymmetric biases > 
5b > 0: when ^ < Ac, only the coexistence state (•) corresponding to a population of A-voters 
and centrists is stable (flows in dashed arrows), while the fixed point associated with centrism 
consensus (o) is unstable. When 5^/(1 + 5^) < ^ < Ac (as in this sketch), there is also an 
unstable coexistence state (o) consisting of B-voters and centrists. When ^ > Ac, the only stable 
fixed point corresponds to the centrism consensus (dotted arrows), (b) Identical biases Sa = 5b = 
5 > 0: the fixed point [a* ,b* ,c* ) (•) associated with the coexistence of the three parties is stable 
when ^ < Ac (flows in dashed arrows) and unstable otherwise (dotted arrows). See text. 



4.1 Mean field dynamics with asymmetric biases 

When 5a> 5b> 0, the mean field equations Q are characterized by three fixed 
points {a,b,c) ^ {a*,b*,c*) with G (1,2,3): 

{alb\,c*,) = (0,0, 1-C), 

{4,bl,cl) = {l-i;-^,0,h, (3) 

OA OA 

{ai,bi4) = {o,i-i:-^,h- 

OB OB 

The non-trivial fixed point (a2,^9,c|) physically exists and is asymptotically sta- 
ble when ^ < Ac- Otherwise, (aj,^j,Cj) is the only stable fixed point when ^ > 
Ac- Furthermore, there is also another coexistence fixed point (03,^3,1:3) when 
< Sb/ {I + 8b), but a stability analysis reveals that it is always unstable (saddle 
point). As illustrated in Fig.llja), the mean field dynamics ^ therefore predicts 
that the system evolves toward the stable fixed point 

This means that in the case of asymmetric persuasion biases with 5a > Sb, when ^ 
is below the critical threshold Ac = 8a/ {I + Sa), the mean field dynamics leads to 
a coexistence state where the most persuasive voters (A's) coexist with centrists. In 
this case, S-voters are absent form the final state. On the other hand, when ^ > Ac 
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the mean field dynamics evolves toward centrism consensus, with a population 
consisting only of susceptible and committed centrists. 



4.2 Mean field dynamics with a identical bias Sa = 5b = S > 

When the persuasion bias toward opinions A and B is identical, i.e. 5^ = 5^ = 5 > 
0, the rate equations (|2]l can be solved exactly. In fact, by adding the equations of 
(|2]l and using a{t) + b{t) = 1 — — c(f), we obtain 

|c(0 = -^-{a{t)+b{t)) = -5(1 - i;-cit)Mt) - (5) 

This equation can be readily solved and by denoting the initial densities (a(0) ,b{0),c 
{x,y,z), when ^ 7^ Ac, one finds: 

^ l-g+/:gexp(g(l-g-g/g)0 
l+7^exp(5(l-C-C/5)0 ' 

where K = ^j^zf ■ When the biases toward A and B are identical, one has ^ ^ f|j ) = 
0, which implies that a{t)/b{t) = x/y is conserved by (|2]i. Using this property to- 
gether with (|6]l, we obtain a{t) = (^7^) [1 - C - c(f)] and b{t) = (_v/x)a(f). Ac- 
cording to the mean field dynamics is therefore characterized by an exponen- 
tially quick approach to the steady state densities {a* ,b* ,c*), where 




< Ac and 5a = 83 = d 



(7) 



and 



a{t) = -b{t)'^a* = -b* = { nSz-Q ^<^' (8) 

y y \o if C > 4o 

as sketched in Fig.[nb). This means that when the persuasion biases are identical, 
and if the density 1^ of centrist zealots is below Ac, the mean field analysis predicts 
the stable coexistence of A and B voters along with centrists. If C > ^c, the final 
state is composed only of (susceptible and committed) centrists. Furthermore, in 
the critical case ^ = Ac, one has the long-time behavior c{t) — (1 + ^ (^0"^ 
when f — > 00. In this critical case one is left with a density = (1 + of 
susceptible centrists along with a fraction 11, of zealot centrists, but no A and B 
voters. 

According to these mean field results, summarized in Fig.[Tl when the fraction 
^ of zealots in the population is low, the dynamics generally reaches a stationary 
state where more than one species coexists. Hence, the mean field analysis predicts 
the transition from a coexistence and "multicultural" phase (1^ < Ac) to a phase 
dominated by centrism (1^ > Ac) at the critical value 1^ = Ac- As discussed in the 
next section, this picture changes drastically when the population is of finite size 
and demographic fluctuations need to be taken into account. 
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5 Stochastic dynamics in a finite population 

When the population is of finite size A'^, demographic (random) fluctuations signif- 
icantly alter the mean field dynamics and ultimately impose the centrist consensus. 
The latter is associated with the system's unique absorbing state that consists of 
N ~ t susceptible centrists (C) and centrist zealots (C^). While the system's fate 
is known, we are here interested in the dependence of the consensus (exit) time 
on the population size and on the parameters 5a, ds and ^. In particular, we study 
the circumstances under which cultural diversity is maintained (long-lived coexis- 
tence), and when consensus is quickly reached. These questions are addressed by 
modeling the evolution with continuous-time birth-death processes. 

In the general case of different persuasion biases (5a > dg), at each interaction 
the number of A and B voters can increase or decrease by one (i.e. j ^ 7 ± 1 and 
^ ^ ^± 1) with rates T'^^^ and T^^, respectively (see ([T|). The population's com- 
position therefore evolves according to a bivariate birth-death process described 
by the probability ^j^kiO that at time f there are j voters of type A and a number 
k of B voters in the population. This probability obeys the master equation 1221 : 



dt 



rf+,,^y-i,,(0 + 7:f;-i,,^,+i,,(f) + 7;^+_i-?3^,M-i(0 



where the state space is bounded since j,k e [0,N — i^] and, to account for the 
fact that j = k = corresponds to the absorbing state with only centrists (C's and 

's), one has T^^ = T^^ = and ^j,k{t) =OfoT j > N -i^ and j < 0, and for 
k> N — and k <0. According to ([T]|, the transition rate T^j^ is associated with 
the reaction of the pair AC^AA and Tf+ with BC BB, while Tf-~ and Tf~ are 
respectively associated with AC /C^ CC/C/^ and BC/C^ — > CC/C^. Since the 
probabilities of picking the pairs AC and AC^ are respectively ji/N{N — 1) and 

/N{N — 1 ), and similarly for the pairs BC and BC^, the transition rates read 

tA+ (, I H J^N-l^-j-k) . KN-j-k) 



N{N-l) ' J^" N{N-l) 

tB+_,. MN-i^-j-k) ,__ k{N-j-k) 



where we have used j + k — N — £^ — i. 

The mathematical treatment greatly simplifies in the case of identical persua- 
sion biases, when Sa = dg = d > 0. In this case, centrists C interact in the same 
manner with A and B voters (see ([T]!) and at each time step their number can 
change by ±1 with transition rates T^, respectively. The dynamics can thus be 
mapped onto a single -variate process in terms of the probability Pi{t) of finding £ 
susceptible centrists in the population at time t obeying the master equation 

^^ = T+,Pe_i{t) + T+,Pe+i{t)-[T+ + Tf]Pfit), (II) 
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where, to account for the absorbing boundary conditions, we impose T^_^^ = 

and Pi{t) = iov I > N — and ^ < 0. According to O, and from the above 
discussion, the interaction XC XX, with X G {A,B}, is associated with the rate 
Tf — T^i^ + T^j^. In the same manner, corresponds to the reactions XY CY, 

with Y € {C,C^}, and thus T^^ = T^^ + T^^ . The transition rates of the master 
equation (1111) therefore read 

T^ = - \ and 7}- = (1+5)——^—^. (12) 

Here we are chiefly interested in the mean consensus time (MCT) T necessary 
to reach the absorbing boundary {j ,k,£) — {0,0,N — £i^) starting from a population 
with initial densities (c((0),/?(0),c(0)) = {x,y,z), while ^ remains fixed. The MCT 
will be studied both for ( 19110b and (II 11121) . In the case of identical persuasion bi- 
ases, the MCT of the birth-death process (111112b will be obtained analytically. The 
direct treatment of the process ( I9I10I) describing the 3CVM with asymmetric per- 
suasion biases is difficult, but its dynamics is exactly reproduced by the Gillespie 
algorithm used in our simulations ||28]| . Furthermore, in Sec. 7 we will use the an- 
alytical results obtained in Sec. 6 for the unbiased case to analytically characterize 
the MCT also when 5a > 5^. 



6 Mean consensus time in tlie 3CVM witli centrist zealots and identical 
persuasion biases 

In this section, we study the mean consensus time (MCT) of the birth-death pro- 
cess defined by ( I11I12I) when the persuasion biases are identical, i.e. 5a = 5b = 
5 > 0, and the population size ^ 1 is large (but finite). We consider the contin- 
uum limit and treat z = £/N and C = as continuous variables. The transition 
rates (fT2l) therefore become 

r+(z) = (l-C-z)(z+?), and r-(z) = (l+5)(l-C-z)z. (13) 
In this setting, the MCT obeys the backward master equation 12211111 

(r+(z) + T-{z)) t(z) =A + T- (z)t(z -A) + T+{z)t{z + A), (14) 

where the time has been rescaled according to t ^ t/N and A = (the time 
increment A hence matches how the density of C's changes at each interaction). 

In principle ( 114b is solvable, but the formal exact solution is an unwieldy and 
non-enlightening expression 0221 . A more tractable and insightful approach is pro- 
vided by system size expansions of (114b in powers of zi . In the realm of the diffu- 
sion theory II27II22I . one Taylor expands (114b to second order in A and obtains the 
backward Fokker-Planck equation (bPPE) 

- T-{z))T'iz) + + T-iz))/2N)T"{z) = -1 
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This bFPE holds for all single-coordinate birth-death processes 1221 and, with 
( fTST l. here explicitly reads 



dz 2N ' dz^ 



= -1. (15) 



This equation is supplemented by the absorbing and reflecting boundary condi- 
tions 'r(l — 1^) = and t'(0) = 0, respectively. The condition t(1 — i^) = ac- 
counts for the fact that the dynamics ends when z = I — ^ (i.e. £ = N — £i^) and the 
centrism consensus (absorbing state) is reached. The reflecting condition t'(0) — 
ensures that the system cannot be "exited" by the end z = (i.e. ^ = 0) ll22l . 
Eq. (fTSl l gives the MCT in the framework of the diffusion theory, which is cer- 
tainly an accurate approximation when the intensity of the deterministic drift is 
much less than the intensity of the fluctuations (oc AT" ^2) |I27'/29'!'301. The formal 
solution of (115b can be obtained by standard methods and reads l,22J : 



t{z)=2N [ dve^'^^'^ f du- 
Jz Jo (l 



Jz Jo (l-C-M)(C + (2 + 5)a)' 

where 



(16) 



^{x) =nx~mln[{2 + 5)x+Q, (17) 

25 ^ 4C(l + 5) 

with and 

The result ( 116b . is valid in both regimes ^ < Ac and ^ > Ac, and is particularly 
relevant when 5 and are of order A^^', and there is an effective competition 
between the deterministic nonlinear drift and the diffusive noise. It is worth em- 
phasizing that in the case of identical biases, centrists do not discriminate between 
A and B voters and the MCT is thus a function ofx + y= 1 — ^ — z (total initial 
density of radical voters) or, equivalently, of the initial density z of C's. A particu- 
larly interesting situation arises when z — and the system consists of a vanishing 
minority of centrists and an overwhelming majority of radical {A and B) voters. 
Below, we show that two very different types of behaviors arise: a form of cultural 
diversity is maintained over a long period of time when ^ < Ac, whereas centrism 
consensus is quickly reached when ^ > Ac. 



6.1 Long-lived coexistence and consensus time 

We first study the situation where < Ac or, equivalently, 5>^/(l — ^).In this 
case, the mean field rate equations (|2]l predict the coexistence of the three opinions 
with a density of susceptible centrists (see Fig.lljb)). Here, by paying special 
attention to the interesting situation where < 5 <C 1 and A'^ ^ 1, we analyze 
how this picture is altered by random fluctuations and how centrism consensus is 
ultimately reached. From the results of stochastic simulations reported in Fig.|2] it 
appears that the MCT decreases monotonically with the density of centrist zealots 
^ and increases with the persuasion bias 5. The MCT also appears to grow steeply 
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Fig. 2 (Color online). Functional dependence of the MCT when ^ < and 5^ = 5/; = 5. 
(a) t(z) as function of the (rescaled) initial concentration of susceptible centrists z/(l — Q. 
Results of stochastic simulations (symbols) averaged over 5 x 10** samples are compared with 
formula ^ (curves) for (5, Q = (0.04,0.02) (x, thick dashed), (0.035,0.02) (A, thin dashed), 
(0.06,0.04) (o, thick solid), and (0.08,0.06) (o, thin solid). The population size is A' = 500. (b) 
Logarithm of the MCT as function oiN with C, = 0.02 and 5 = 0.04(o) and 5 = 0.03(A) starting 
from a system without susceptible centrists: log^ t(0) has been computed using stochastic sim- 
ulations (symbols, averaged over 1000 samples with N = 100 — 2000), and is compared with the 
pred ictions of l fT6t (solid) and the WKB results J26) (x). The (rescaled) asymptotic predictions 
( l20t are shown as dashed lines. 



(almost exponentially) with the population size A'^, see Fig.|2jb). Moreover, the re- 
sults of Fig. |2ta) indicate that T t(0) when the initial number of susceptible 
centrists is small (z ^ \ — 11,), which implies that in such a regime the MCT is 
essentially independent of the initial density z. Below, these observations are sub- 
stantiated by the analysis of the functional dependence of the MCT on 5, and A'^, 
first using the Fokker-Planck equation and then with the WKB method. 



6.7.7 Mean consensus time with the Fokker-Planck equation when 11, < Ac 

In this section, we analyze the predictions of ( fT6] l when 11, < and z ^ 1, and 
obtain the leading contribution of the mean consensus time in the framework of 
the Fokker-Planck equation ( 115b using a saddle -point approximation. 

When A'^ ^ 1 and ^ < Ac, the function exp{—N^{u)) has an isolated peak at 
u* = 1^/5 and one can Taylor-expand ^(m) around u* which yields exp(— A^^(m)) = 
exp(— TV [^{u*) + (m — m*)^^"(m*)/2] . We can use this expansion to evaluate the 
inner integral of the right-hand-side (RHS) of ( fT6] l with a saddle-point approxima- 
tion, yielding 



lo ""(l-C-")(C + (2 + 5)« 



, ^ r du e- W2)(«-?/^)'-^"(?/^) . 

J — oo 



In setting the boundaries to ±°o we have assumed that the peak is sufficiently 
separated from the absorbing boundary, i.e. we henceforth assume that C/^ is 
well separated from the values and 1 . By performing the Gaussian integral, the 
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expression (fT6] l becomes 



As ^(v) is an increasing function ofvon i^/5<v< 1 — the main contri- 
bution to the integral on the RHS of ( 1181) arises from v ~ 1 — and therefore 
ri-Crf^g/v,ir(v) ^ j.i-?^^^iv[,.^(i-f)-(i-c-.),^'(i-?)] _ eiv,^{i-?)/^. Hence, with 
(I18l l. and provided that the initial population does not mainly consist of suscepti- 
ble centrists (i.e. if z ^ 1), the leading contribution to the MCT in the realm of the 
diffusion approximation is 

T ^ TV-i/ZgWrls,?) ^^-i/2gyv(,i?(i-?)-^(?/5))^ (19) 

where the functional dependence of the exponent y(5,i^) = ^(l — — 
is illustrated in Fig.[3j 

As we are particularly interested in the limit of weak persuasion biases and 
small fraction of zealots, i.e. ^ < 5 <C; 1 (with z ^ 1 — 0' th^ expansion of the 
exponent /yields y( 5, = 5- ^(1+ In (5/0) + i^(5^). With (19^, one therefore 
obtains the following asymptotic result when ^ 1 : 

T^7v->/2.^(^-?)(ly^ (20) 

The concise result ( 120b provides the leading contribution to the MCT in the realm 
of the diffusion theory and is valid as long as < 5 ^ N~^l^ . According to this 
result, the leading contribution to the MCT grows exponentially with the popu- 
lation size and with the difference 5 — (in a manner that is independent of the 
initial condition z, see below), and decreases with The exponential growth of 
the MCT with N{d — ^) is confirmed by the stochastic simulations reported in 
Fig. |2l where we find an excellent agreement with the predictions of ( 116b . The 
results of Fig. [2a) also show that the MCT depends "weakly" on the initial den- 
sity z of susceptible centrists: the MCT significantly deviates from the value t(0) 
only when z/(l — Q > 0.7, i.e. in the special case where there is already an over- 
whelming initial majority of centrists. In Figure |2b), we see that the asymptotic 
behavior of the MCT when N{5 ~ Q I is aptly captured by ( l20l) . 

Furthermore, to understand the influence of the term (O 5)'^^ on the RHS of 
( |20] |. it is worth considering the transformation (5, — > (5', 0) = (5 + a, ^ + a), 
where a is a given (small) real number. Under such a transformation, the exponen- 
tial term of ^ is left unaltered, but {^Sf^ ((C + a)/(5 + a)f^^+"^ and 
therefore the MCT decreases if a > and increases otherwise. This nontrivial in- 
terplay between 5 and ^ is confirmed by the results of Fig. |2f a) and is illustrated 
in Fig. [3] where we see that the exponent / decays along the line d — ^ =constant. 

6.7.2 Metastability and WKB treatment when N{5 — 0^1 



The asymptotic behavior ( l20b of the MCT has been derived from the solution of 
the Fokker-Planck equation ( 115b . by assuming that < 5 ^ 1 . In this section, we 
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according to il9\ : 7 grows and decreases monotonically with 5 and ^, respectively. 



use a complementary approach in terms of the Wentzel-Kramers-Brillouin (WKB) 
theory 131113211 to obtain the MCT to the next-to-leading order in the limit where 
N{5 - p > 1 when ^ < A^. 

While the mean field analysis Q predicts the stable coexistence of all opinions 
in this regime, when the population is large but finite, the system lingers around the 
coexistence state (with a density of C's) before centrism consensus is eventu- 
ally reached. The dynamics is thus characterized by metastability, a phenomenon 
that is well described by the WKB method f3P|. The WKB approximation is an 
asymptotic theory frequently used in the semi-classical treatment of quantum me- 
chanics Bill , and t has also been used to study stochastic processes |[32l . Recently, 
there have been numerous applications of the WKB approach to study problems 
of metastability arising in population dynamics, as well as in evolutionary games 
and population genetics, see e.g. |29 33 30 1. As the WKB method has never been 
used in the framework of opinion dynamics to the best of our knowledge, we here 
follow Refs. II29II301 and outline the essence of the method and show how it can 
be used to compute the MCT when N{8 — 1. 

The WKB treatment is based on a size expansion of the master equation 
dllb using an exponential ansatz for the probability distribution (see ( I23l l and 
Ref. f3ri ). It relies on two main assumptions: (i) the population size is large (i.e. 

I) and one can work in the continuum limit; (ii) the system quickly relaxes 
toward the metastable state and, from there, reaches the absorbing (consensus) 
state after a time growing dramatically with A'^ Il32ll33l|29l l30 |. Here, this implicitly 
means that we assume that the initial density of centrists is low enough (z <C 1 — Q, 
to ensure that the metastable state is always reached before consensus is attained 
(see Fig.[Tli. In this setting, the mean consensus time Twkb is the decay time of the 
metastable state described by the centrism quasi-stationary probability distribu- 
tion (QSD) Tif. The QSD is obtained by approximating Po<(<fj^(^ 7r(?e~'/^wKB 

and Pn-c^ — 1 — e^'/^wB ||29II30I . Using these expressions in the master equation 
(lllb . one obtains the quasi-stationary master equation 

T+^ne-i + Tf^^ne+i - [T+ + Tf] m = 0, (21) 
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where an exponentially small term tt^/twkb has been neglected and we have used 
the transition rates (112b . The decay time of the metastable state coincides with the 
MCT and is given by the flux of probability into the absorbing state, i.e. 



dT+{l-Q 



dz 



N ' 



(22) 



where we have used the continuum limit z = i/N, = Knz, = 7z{z) and 7]^_£^_j = 

|iir+(l - 1^) /dz\/N, with the continuous rates ^ II29II30II . According to ([23, the 
calculation of T requires to solve the equation ( 1211) for the QSD. This can fruitfully 
be done using the WKB approach 132], which is based on the ansatz 

7r(z)^^e-'^^W-^iW, (23) 

where S{z) and (z) are respectively the system's "action" and "amplitude", 
while jz/ is a normalization constant. The action is found by substituting i23i 
into (1211 and by keeping the leading order in A^, which yields T'^{z)[e^ — 1] + 
T- (z) [e-^'(^') - 1] = 0, with S' = dS/dz Il32ll33ll29ll301. With ([I3]l, one thus finds 

Siu) = - l"\n[T+iv)/T-{v)]dv = uln (^^^^^ ) - Cln(«+ 0- (24) 

The constant £/ is determined by normalization of the Gaussian approximation 
7i{z) - ^e-'VS(z*)-(A'/2)5"(z*)(z-z*)' around the metastable state z* = C Z^- A Gaus- 
sian integration yields ^ ^j^j^ leading order, n{z) e~'^('^(^)~'^(?/^)). 
With ( 122124b . the leading contribution to the MCT therefore reads O3ll29|[30l : 

It is worth noting that, to low orders in d and ^, the exponent Ywkb of ( 125 b reads 
1^vkb(5,C) ^ 5(1 - Q-S{U5) = (1 - C)ln(l +5)(1 - Q + ClnC(l + 

= d C(l +ln(5/C)) - - C2)/2 + ^(5^). 

Hence, the exponents y of the Fokker-Planck result ( fTOl l and ]^vkb of the WKB 
treatment ( |25] | coincide to first (leading) order in < 5 <C 1, but their next-to- 
leading-order contributions are different. 

The calculation of the subleading correction to the MCT ( 125) is feasible but 
more involved. In fact, it necessitates (z) and requires to match the ansatz ( 1231) 
with the solution of the (121b linearized about the absorbing boundary z = 1 — C- 
Here, we quote the final result and refer to Ref. 12911 for a detailed calculation and 
a related discussion 0: 



N / 9- \Ni; 



(25) 




Twkb:^W-(^1 ^s{l-Q-Cr \5U~T)) ■ ^^^^ 

One slight difference with the generic result given in Ref. 1 29 1 lies in the fact that here the 
absorbing boundary is at z = 1 — ^. Furthermore, in the final result l l26t we have chosen the 
same timescale as in Sec. 6.1.1 and divided <22t by A'. This allows a direct comparison with the 
results of the Fokker-Planck treatment and differs by a factor A'^' from that of I29I30I . 
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Fig. 4 (Color online). Functional dependence of the MCT when [, > Ac and 5^ = 5^ = 5. 
(a) t(z) as function of the (rescaled) initial concentration of susceptible centrists z/i l — Q for 
various values of 5 and ^. The results of stochastic simulations (symbols), averaged over 5 x 10* 
samples are compared with formula <16t (curves) for {5,Q = (0.04,0.06)(x, thick dashed), 
(0.06,0.08)(o, thick solid), (0.08,0.10)(o, thin dashed) and (0.035,0.08) (V, thin solid), (b) The 
MCT as function of logjg A', with A' = 100 — 10000: starting from a system without susceptible 
centrists with (5,Q = (0.06,0.08) (o) and (5,Q = (0.08,0.10) (o), -r(O) has been computed 
using stochastic simulations (symbols) and the formula U6t (solid/dashed lines). 



As illustrated by Fig.|2tb), when the population size is sufficiently large and N{d — 
Q ^ 1, ( l26l) is in good agreement with the predictions of the FPE and with the 
results of stochastic simulations. Two additional remarks are in order: 

i. In contrast to the prediction ( 116b of the diffusion approximation, the results 
( 125126b for Twkb only depend on the parameters d and ^ but are independent 
of the initial density of centrists z, and are thus to be compared with t(0). 
This is because the MCT is here computed /ro/w the metastable state, which is 
certainly always legitimate when the initial population does not consist of an 
overwhelming majority of centrists (i.e. provided that z, ^ 1). 

ii. The accuracy of the WKB approximation improves when A^(5'(l — i^) — 5(1^/5))^ 
1, see e.g. II29II30I . as illustrated in Fig.|2lb) where it is shown that the WKB 
results are more accurate when 5 = 0.04 (and A'^ > 1000) than for 8 = 0.03. 



6.2 Centrism consensus time when ^ > Ac 

When 1^ > Ac, the rate equation (|6| is not characterized by an interior fixed point 
and the mean field rate equations ([2]l predict that centrism consensus is quickly 
reached (after a short coexistence of the three opinions). This mean field scenario 
essentially also holds in a finite population when > as illustrated by Fig.|4la), 
where at fixed it is found that 

- the centrism consensus time T grows when the density z of susceptible centrists 
decays; 

- T decreases monotonically when 5 decreases and increases; 

- at fixed value of ^ — 5, T increases monotonically when 5 and decrease; 

- the parameters and 5 appear to have only a weak effect on the MCT when 

z^O; 
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- the predictions ( fT6l) obtained from the backward Fokker-Planck equations are 
in excellent agreement with the results of stochastic simulations. 

When ^ > Ac, the integrals in ( 116b are difficult to evaluate analytically since the 
function sxp{—NJ^{u)) has no isolated peaks and thus no saddle-point approxi- 
mation can be carried out. Nevertheless ( 1161) can be studied numerically and it is 
found that when z — !■ 0, i.e. when the initial population consists only of A and B 
voters (along with centrist zealots, but no susceptible centrists), one has 



This result, compatible with the above points (iv) and (v), is confirmed by Fig.|4fb) 
where, when is large and increases, it is shown that 

- T grows linearly with logjgA'^, with a slope that appears to slightly increase 
when ^ and d decrease, as shown in Fig.lHb) when > 1000; 

- the formula ( 116b fully captures the dependence of T on A'^ for all population 
sizes in the range A' = 10^ — 10"*. 

7 Mean consensus time in tlie 3CVM witli centrist zealots and asymmetric 
persuasion biases 

When the biases toward A and B voters are different, i.e. when 5a > 5b > 0, the 
system's dynamics is described by a bivariate birth-death process characterized by 
the master equation (|9]l. The mean time to reach centrism consensus starting with 
a population comprising a fraction of x voters of type A and a density y of S-voters 
thus obeys the following backward master equation II22I : 



z{x,y)=A + T-T{x-A,y) + T+Tix + A,y) + T-T{x,y-A) + T+T{x,y + A) 



where A = N . Since we consider a a finite but large population size (A^ ^ 1), 
the transition rates T^y are the continuous versions of dlOl) . i.e. T-^ = (1 + 5a)xz, 
T~ = jc(z+ 1^), 7"^+ = (1 + 5B)yz and T~ = y(z+ 1^). It is generally difficult to solve 
(128b and even a size expansion is not very useful to make analytical progress. Yet, 
the bivariate master equation (128b can be simplified in the limits where x — > 
and y — > 0, and the asymptotic behavior of the MCT can be obtained when the 
population size is large (see below). In fact, when y = and x = 1 — i^— z>0, 
= while = (z) and therefore t{x, 0) = Tg^ (z = 1 — ^ —x), where Tg^ 
is the solution of the single-variate master equation (1141) . with 5 replaced by 5a. 
Similarly, when x = and y=l~^— z>0, one has T(0,y) = (z = I — ^ —y). 
Building on our knowledge of the case 5a = 5b (Sec. 6), we know that the solution 
of (114b increases monotonically with 5 when ^ is kept fixed (see Fig. |2]l. Since 
5a > 5b, this implies t(x, 0) = Tg^ (z = 1 - C - Jf) > T:{0,y) = T^^ (z = 1 - C - .v)- 
More generally, at fixed values of z and ^, the largest MCT is attained when the 
population of non-centrists consists of a maximum number of most persuasive 
radical voters, i.e. N — I — £^ voters of species A. Similarly, the shortest MCT at 
fixed values of z and ^ is obtained when the initial sub-population of non-centrists 
only consists of N — £ — l/^ voters of species B. Hence, the MCT respectively 



T ^ InN, when A^ » 1 . 



(27) 




(28) 
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0.2 0.4 0.6 0.8 500 looo isoo 2000 

-^/(1"C) with xjy fixed ^ 



Fig. 5 {Color online). Functional dependence of the mean consensus time when if < 4^ and 
5a > with f = 0.02. (a) The MCT as function of the rescaled initial density z/(l — if ) with 
the ratio xjy held fixed, see text. Here, = (0.04,0.03) withx = y (o) w\Ay = 1.5x (x), 

the population size is N = 500. The results of stochastic simulations (symbols) averaged over 
2.5 X 10* samples are compared with the results obtained for 5a = Sg = 0.04 (A, so lid) and 
Sa = Sb = 0.03 (o, dashed). The analytical results (solid and dashed lines) are based on l l30l31t . 
(b) Logarithm of the MCT as function of A' (with N = 100 — 2000) for the same parameters as 
in (a): logjoT(x + y= 1 — f ) computed usin g stoch astic simulations (o, x) averaged over 10^ 
samples is compared with the predictions of ( 13013 U for 5 = 0.04 (solid) and 5 = 0.03 (dashed). 



increases and decreases with the initial fraction of A-voters and S-voters in the 
population, and the quantities Tg^ and Tg^ are upper and lower bounds of the MCT 
[see Eas (l30l3TT) .1: 

Z5,il-^-x-y)< zix,y) < 15^(1 -^-x-y), (29) 

as illustrated by Figs. |5] and |6] see below. According to this discussion, we no- 
tice that a quantitative difference between the cases of asymmetric and identical 
persuasion biases lies in the influence of the initial condition: when Sa 7^ 5b, the 
MCT depends on the separate initial densities x of A's and y of S's (not only on 
their sum, as in the symmetric case with Sa = dg)- In what follows, we show that 
the asymptotic behavior of the MCT is independent of the initial condition when 
N{dA — 8b) ^ 1 and significant insight into its properties can be gained from the 
results of Sec. 6 combined with the mean field analysis of Sec. 4. 



7.1 Long-lived coexistence and consensus time when 1^ < Ac and 5a > 5b 

When ^ < Ac = 5a/ {I + 5a), the mean field analysis (Sec. 4) has revealed the 
asymptotic stability of the state where A-voters and centrists coexist. As seen 
in Sec. 6.1, when the population size is large but finite, we expect the coexis- 
tence state to be long-lived and to decay after a time growing exponentially with 
the population size A'^ (to leading order). If the population initial composition is 
{x,y,z) = (1 — 1^,0,0) and only consists of A-voters and centrist zealots, one re- 
covers a single-variate process (as in Sec. 6) and the MCT is z{x = 1 — i^,0) = 
Tg^iz = 0), where, in the realm of the Fokker- Planck equation, Tg^ (z) is directly 



21 



inferred from (fT6l) using 

Tsiz) =2N dv .^^^(W / du—- . ^. (30) 

Jz Jo [l - Q -u){Q + {2 + d)u) 

with 

Similarly, if the initial state is (ji:,y,z) = (0, 1 — CjO), one has t(0, 1 — i^) = ^^^(O). 
Since Sa > 5b and T5 - Af- VZgA'rls.?), with an exponent 7= ^^(1 - C) - ^5(C/5) 
that grows monotonically with 5 when ^ is fixed (see Sec. 6.1 and Fig.O, one has 
t(x,0) = T5^(0) -Af-VZgA'rlSA,?) > t:(o,3;) = Tg^[0) when 7V(5a - 5b) > 1 and 
the initial density of susceptible centrists is vanishingly small (z — ^ 0). 

In the case of a population with an arbitrary initial density {x,y,z) of A,B 
and C's, the system almost surely reaches the metastable state (c!2,^2''-2) ~ 
^ — 1^/54,0,1^/54), as prescribed by the mean field equations (|2]l [see Q and 
Fig.lifa)], before reaching centrism consensus. Hence, the less persuasive of the 
radical voters (B's) quickly disappear from the population and the metastable state 
only consists of centrists and voters of the most persuasive party (A's). Apart 
from the short transient necessary to reach the metastable state, the MCT is thus 
obtained from the mean time to reach consensus starting from the fixed point 
((221 ^2' "-D where the population comprises only A-voters and centrists. In this 
situation, when N{8a~ ds) ^ I, the leading contribution to the mean consensus 
time T(x,y) is obtained from a single-variate backward master equation [like ( 114b 
of Sec. 6] and is given by Tg^ (z) with ( 1301) . In particular, when the initial density z 
of susceptible centrists is very small (i.e. z ^ 1 — 1^) and A'(5a — 5b) ^ 1, one has 
T{x+y ~ 1 — 1^) ~ Tg^ (0) and, from the results ( l20l) and ( l25T l, we infer the leading 
contribution to the MCT when 5b < 5^ ^ 1 and 1^ ^ 1 : 

/ r \ Ni 

r{x + y^l-Qc^Ts,iO)^N-'/^e^^'^-'^^^-l-j . (32) 

The results of stochastic simulations reported in Fig.|5]corroborate the above anal- 
ysis. In fact, Fig.|5ta) reveals the dependence of the MCT on the initial conditions: 
T grows monotonically with z^^ and x/y, but exhibits only a weak dependence on 
x/y. Moreover, the comparison with the case of identical biases illustrates that 
Tsg{z= I — ^ —x — y) < T{x,y) < Tg^ (z = I — ^ —x — y) which confirms ( |29] |. 
Fig. |5fb) shows that asymptotically the MCT grows exponentially (to leading or- 
der) with the population sze when Af(54 — 5b) ^ 1 and z — 0, in full agreement 
with the behavior z{x+y = 1 - ^) t^^ (0) predicted by (|20]l-(I22]l. The figure|5];b) 
also illustrates that the asymptotic behavior of the MCT is independent of the ratio 
x/y and therefore of the detailed composition of the initial population. 

7.2 Mean consensus time when > Ac and 5a > 5b 

When ^ > Ac, the only stable fixed point corresponds to centrism consensus. 
From the analysis of Sec. 6.2, we know that in a finite population with initial 
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z/{l-Q with x/y iked logio 

Fig. 6 (Color online). Functional dependence of the MCT when ^ > A^. and 5a > 5^, with 
^ = 0.08. (a) The MCT as function of the rescaled initial density z/{l — Q with the ratio x/y held 
fixed, see text. Here, {5^,5^) = (0.04,0.03) with.i: = >> (o) and y= 1.5x (x), the population size 
is N = 1000. The results of stochastic simulations (symbols), averaged over 3.5 x lO'* samples, 
are compared with results obtained for Sa = 5b = 5 = 0.04 ( A, solid) and 5 = 0.03 (o, dashed). 
The analytical results (solid and dashed) are based on l l30l31t . (b) T{x + y= 1 — if ) as function of 
Nina semi-logarithmic scale for the same parameters as in (a). Results of stochastic simulations 
(symbols) for A' = 100 — 10000 are compared w ith Tg„ (0) (thin solid) and ig^ (0) (thick and, as 
a guide for the eyes, dashed line) obtained from <30l31t . see text. 



composition {x = 1 - ^,0,0), the MCT is t(1 - i^,0) = Tg^iz = 0) and can be 
obtained from ( I30I) . Similarly, for an initial composition (0,y, 0), the MCT is 
t(0, 1 — i^) = Tgg (z = 0). When N ^ 1 and the initial number of susceptible cen- 
trists is small, i.e. z=l — i^— jc — y— >0, we have seen in Sec. 5.2 (see Fig.|4tb)) 
that Tg^ (0) ^ Tgg (0) ^ InA'^. Furthermore, according to ( 1291 ) Zg^ and T^^ are upper 
and lower bounds of the MCT. The results of Sec. 6.2 hence imply 

(0) < t(jc + y ~ 1 - 1^) < Tg^ (0) , and therefore 
T(x + }'~l-C)'-T5^(0)-T5a(0)-lnAr, whenA^>l. (33) 

These findings are corroborated by the results of stochastic simulations re- 
ported in Fig. |6] that confirms ( |29] |: Tggiz = I - -x-y) < T{x,y) < Tg^iz = 
I — ^ —x — y). The dependence of the MCT on the initial conditions is illustrated 
in Fig. 13 a) where T{x,y) grows monotonically with z~^ and x/y (held constant). 
We notice the dependence on x/y is particularly weak. Fig.|6tb) confirms the log- 
arithmic growth with A'^ of the MCT when the population size becomes very large, 
in full agreement ( 133b . In particular, the symbols in Fig. [S^b) are found to be 
aligned parallel to the thick and dashed lines when N ^ 1, which implies that 

asymptotically e^(-'^^>') Q{x,y)e'^^A^^^ or T{x,y) '^5^{0) +ln2(jt,y), where the 
function Q depends weakly on the initial densities {x,y). 



8 Conclusion and discussion 

After having reviewed the main properties of the paradigmatic voter model and of 
some of its generalizations accounting for zealotry, commitment and incompatibil- 
ity (Section 2), we have studied how a committed minority may resist a persuasive 
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majority and how the resuhing competition influences the maintenance of cultural 
diversity. Here, this line of investigation has been carried out in the framework 
of the three-party constrained voter model with "centrist zealots". In fact, follow- 
ing the basic idea of the bounded-compromise and three-party constrained voter 
models, we have considered that individuals of two radical parties (A and B) are 
incompatible among them and compete to impose their consensus by convincing 
centrists that, in turn, can convince A's and B's. To reflect the fact that radical (A 
and B) voters are more persuasive than centrists, we have considered that there 
are persuasion biases 5a and 5b toward opinions A and B (with 5a > 5b)- Fur- 
thermore, we have also assumed that the population includes susceptible centrists 
(C) and a small fraction ^ of centrist zealots (C^). Whereas C's can radicalize and 
adopt either the state A or B, centrist zealots are committed individuals that always 
remain centrists and oppose any consensus that is not centrism. Hence, the persua- 
sion of A's and B's is resisted and opposed by the commitment of centrist zealots. 
This results in a subtle competition between commitment and persuasion that has 
been studied in the mean field limit and in a well-mixed population of finite size 
A' (on a complete graph). 

In an infinitely large population, there is a continuous transition between a 
coexistence ("multicultural") phase, that is stable when the fraction of centrist 
zealots is below the critical threshold Ac = 5a/ {I + 5a), and a phase where cen- 
trism prevails when ^ > A^. This scenario changes in finite populations when 
demographic fluctuations ultimately lead to centrism consensus. The competition 
between commitment and persuasion is thus characterized by the mean consensus 
time and has been investigated in terms of single- variate and bivariate birth-death 
processes. The case of identical persuasion bias has been thoroughly analyzed and 
the mean consensus time has been computed in the realm of the diffusion theory 
(Fokker-Planck equation) and with the WKB method. In a large population, it 
has been shown that the long-time dynamics of the case with asymmetric biases 
(Sa> Sb > 0) can be described in terms of a single- variate birth-death process 
involving only centrists and A-voters. Hence, the asymptotic behavior in a large 
population of the mean consensus time and its upper and lower bounds have been 
obtained analytically also when 5a > 5b- In fact, the main difference between 
the cases 5a = 5b and 5a > 5b lies in the composition of the metastable state: it 
is characterized by the coexistence of A-voters and centrists (no B-voters) when 
5a> 5b, whereas voters of the three parties coexist when the persuasion biases 
are identical {5a = 5b)- 

Our findings, can be summarized as follows; When the fraction of centrist 
zealots in the population of size is low (^ < Ac), we have shown that the in- 
terplay between commitment and persuasion restdts in a metastable coexistence 
(multicultural) phase. The mean consensus time thus corresponds to the decay rate 
of the metastable state and, when centrists are initially in the minority, it has been 
shown to grow as T ~ Af~'/^e'^^(^'*'^) with an exponent 7 determined explicitly. 
When the biases are asymmetric, the less persuasive of the radical voters (B's) 
quickly disappear from the population and the metastable state only consists of 
centrists and voters of the most persuasive party (A's), whereas opiiuons of all 
parties coexist when 5a = 5b. On the other hand, when C, > Af., the system attains 
centrism consensus in a time that scales logarithmically with the population size, 
i.e. T ~ InA^. These features are found to be robust since they neither depend of the 



24 



detailed initial composition of the population nor on the weakest persuasion bias. 
It is worth noticing that a similar type of behavior concerning the mean fixation 
time has recently been found in evolutionary games, see e.g. Refs. I30II34II . 

In conclusion, while the mean consensus time of the three-party constrained 
voter model without zealots scales linearly with A'^, we have here shown that the 
presence of small fraction of zealots and the existence of persuasion biases results 
in a rich dynamics characterized either by a prolonged maintenance of a multi- 
cultural phase, or by a quick realization of centrism consensus. The long-lived 
coexistence state generally consists of voters of the party supported by the most 
committed individuals along with the most persuasive voters. Since the descrip- 
tion of the (meta-)stability of the multicultural phase and the mean time to reach 
the consensus cannot be obtained solely from the mean field rate equations, this 
work illustrates the relevance of statistical physics methods to study problems of 
"opinion dynamics". 
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